Constraints on the Quasiparticle Density of States in High-T c Superconductors 
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In this Letter we present new tunneling data on YBa2Cu307 thin films by low tempera- 
ture scanning tunneling spectroscopy. Unusual peak-dip-hump features, previously reported in 
Bi2Sr2CaCu208+i , are also found in YBa2Cu307. To analyse these common signatures we pro- 
pose a new heuristic model in which, in addition to the d-wave symmetry, the gap function is energy 
dependent. A simple expression for the quasiparticle density of states is derived, giving an excellent 
agreement with the experiment. The dynamics of the quasiparticle states and the energy scales 
involved in the superconducting transition are discussed. 



Getting to the heart of the quasiparticle density of 
states (DOS) in cuprate superconductors is an active goal 
of current theoretical and experimental research. Indeed, 
in the case of conventional materials, the quasiparticle 
DOS in the superconducting state contains the key infor- 
mation on the pairing. For this reason various methods, 
such as angle-resolved photoemission (ARPES), tunnel- 
ing (TS) and Raman spectroscopies, as well as optical 
conductivity, are extensively used to elucidate the pairing 
mechanism in cuprates. Among these, scanning tunnel- 
ing spectroscopy (STS), with its high spatial and energy 
resolution, emerges as the technique of choice. 

Using such a local probe at low temperature, repro- 
ducible vacuum tunneling spectra have been obtained 
on many cuprates pj-H|, in particular P^S^CaC^Os+a 
(BSCCO) and YBa 2 Cu 3 7 (YBCO). However, one still 
cannot interpret the experimental curves precisely. The 
shape of the DOS in BSCCO has a non-trivial depen- 
dence on doping ff7p||To|] , temperature J7p|-]T6[] and mag- 
netic field |17|Jiq| . The state density is characterized by 
its large value at the Fermi level, unusual quasiparticle 
peaks and dip-hump structures at higher energies (see 
Fig.l). While the T c vs doping curve is bell-shaped, 
the gap width decreases linearly, as the doping increases 
from underdoped to overdoped regimes, with no signif- 
icant change in the DOS shape py£|jic||. The tempera- 
ture dependence is also unconventional: The gap does 
not vanish at the critical temperature and a pseudogap 
persists at T > T c Finally, a low tempera- 

ture pseudogap remains when the superconducting order 
is destroyed, within the vortex core Q , or due to disor- 
der All these features are difficult to explain with a 
BCS-type theory. 

In this Letter we report our new low-temperature STS 
data on YBCO thin films. We observe strongly pro- 
nounced peak-dip-hump structures and a significant state 
density at zero bias. With previous data on BSCCO, our 
observations suggest a common mechanism for these fea- 
tures, linked to the superconducting state in cuprates. In 
order to clarify this effect, we propose a new quasiparticle 
DOS based on an energy-dependent gap function. A par- 
ticular form for this function is inferred from the data and 



has, in addition to <i-wave symmetry, a single minimum 
at a characteristic energy near the gap value. Using such 
a simple approach, all features (wide quasiparticle peaks, 
dips, humps and zero-bias) in the spectra of both YBCO 
and BSCCO may be properly fitted with remarkably few 
parameters. This procedure then allows one to follow the 
quasiparticle states dynamics in the phase transition, in- 
cluding the pseudogap. Moreover, the scaling of the dip 
position with the gap value at different dopings [Q,^ 10 
is a natural consequence of the approach. 
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Fig. 1. dI/dV(V) spectra on BSCCO taken at T = 4.2 K 
(dots). The positions of the quasiparticle peak, dip, hump 
and zero bias, are indicated by P, D, H and ZB, respectively. 
Solid line: the fit by our heuristic DOS (Ao = 46 meV, <5o = 17 
meV, Ao = 0.45, T = 2.9 meV). 

Inset a: Experimental spectrum (dots) and d-wave fit (solid 
line) with an anisotropic normal Fermi surface as in Bot- 
tom curve: difference between experimental and d-wave fit. 
Inset b: Superconducting DOS (solid line) and low tempera- 
ture pseudogap (dotted line), measured in disordered BSCCO 
at T = 4.2 K. 

It is now clear that, in a first approximation, the case 
of BSCCO is consistent with a <i-wave gap function (in- 
set (a) in Fig.l). However the dips and humps, observed 
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beyond the gap in the quasiparticle DOS (Fig.l), cannot 
be simply explained within this approximation. These 
features, existing at low temperatures, disappear at T c 
j7|,PUlO|] , and thus should be related to the supercon- 
ducting state. Using STS, Renner et al. || report on 
normal tip-superconductor (SIN) junctions, and find the 
dip to scale with the gap width for different doping lev- 
els, conserving its characteristic energy near — 2A. This 
behaviour was confirmed in recent reports [0,[lO| where 
both superconductor-vacuum-superconductor (SIS) and 
SIN junctions with BSCCO were studied, and for a large 
range of carrier concentration. In the latter work, the 
dips were found at both negative and positive biases and 
at the energy near ±3A for symmetric SIS junctions and 
at ±2A for SIN ones. Similar behaviour was observed in 
ARPES measurements 0121. 
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Fig. 2. Normalized tunneling spectra of YBCO at T = 4.2 
K (dots). Solid lines: the fit by our new DOS function 
(A = 37 meV, 8 = 21 meV, A = 0.5, T = 5.5 meV for 
spectrum 1; Ao = 26 meV, Sl = 14 meV, Sr = 11 meV, 
A L = 0.7, Ar = 0.45, r = 8 meV for spectrum 2). 
Inset: Raw tunneling spectrum (left curve) observed in the 
same thin film of YBCO at T = 4.2 K. Right curve: The 
same spectrum, normalized to the parabolic background. 

The quasiparticle DOS in YBCO was found to be dif- 
ferent. Compared to BSCCO, a higher density of states 
at the Fermi level (30-60% zero-bias value), a significant 
parabolic background and double peaks, were observed 
Jl7|] . However, no clear evidence for dip signatures were 
reported. We indeed observe spectra of this type in some 
locations of our thin films |r| (inset in Fig. 2). However 
in most of the sample surface a second type of spectrum 
is found (Fig. 2). The doubled peaks are absent there and 
a single broad peak is followed by strongly pronounced 
dips at both occupied and empty states. Another feature 
is an occasional anomaly at roughly +5 meV, similar to 



the resonant state observed in BSCCO and attributed to 
a point defect |2^,|l|. Comparing experimental curves, 
the shapes of our spectra in YBCO p2j and BSCCO are 
quite similar. These data clearly show that the large 
peaks, the dip-hump features and the important zero- 
bias background, are three signatures of the supercon- 
ducting state in cuprates, and are not just a peculiarity 
of BSCCO ]23| ]. None of these features are explained 
within a BCS-type weak coupling approach, even with 
some particular symmetry of the gap function. We now 
describe a new quasiparticle DOS, fitting in detail the 
tunneling spectra, and its consequences. 

Consider a dispersion law of the general form: = 
\A k + A£, where e k is the normal state spectrum. In 
the case of a mean-field superconductor, the gap Ak is 
either constant (strict BCS) or k dependent (BCS-type). 
In the latter case, of particular interest is the <i-wave 
form: A(p) ~ Aocos(2(^), tp being the angle in the a-b 
momentum plane. Then, in polar coordinates, the super- 
conducting tunneling DOS, N S (E), is written in terms of 
a spectral density n s (E, tp): 



N a (E) 



n s (E, ip) dip 



(1) 



Neglecting the angular dependence of the Fermi sur- 
face, the detailed conservation of states n n (e) de — 
n s (E, tp) dE leads to the expression: 



n s {E,p) = n n (E F ) 



E 



y/E*-A(<p)*' 



(2) 



where n n (e) « n n (Ep) = N n (Ep)/2Tr is assumed for 
the normal DOS. Accounting for the anisotropy of the 
Fermi surface leads to the small correction: n n (Ep , p) 
P). Finally, including some pair- breaking is equivalent 
to replacing E — > E — iT and then taking the real part. 
As is well known, n s (E, p) is singular at E = A(tp) and, 
when integrated over tp, gives the famous d-wave DOS 
widely used in the interpretation of the tunneling data. 
Having no other singularities, expressions of this general 
type cannot describe the experiment. 

We go a step further by now assuming that the gap 
function depends on the quasiparticle states, while the 
dispersion law remains unchanged. A simple choice is of 
the form: A(p) — * A(Ek,p), i.e. which ignores any ex- 
plicit momentum dependence. The detailed conservation 
of the number of states follows analogously, giving a new 
expression for the superconducting spectral density: 



n s (E,p) = ^-N n {E F ) E 
Ztt 



A(E,p) 



dA(E, v ) 
dE 



y/E?-A(E,<py 



(3) 



Upon integration (1), the first term of the DOS is con- 
ventional, while the second term has a new and direct 
dependence on the gap function. This result should not 
be confused with strong-coupling, where the second term 
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in (3) is absent. In fact, this formula has been known for 
some time and was derived by J.R. Schrieffer |24| ]. It does 
not apply to conventional superconductors, in which re- 
tardation effects are crucial, and was therefore dropped. 

We now focus on aspects of the spectra which deter- 
mine the precise form of A(E, tp). Recent experimen- 
tal [7| JTl|Jl^| and theoretical |^5| works have focused on 
the dip- hump features in BSCCO. However, we stress 
that the tunneling spectra of both BSCCO and YBCO 
are also characterized by unusually large quasiparticle 
peaks (much larger than the energy resolution of STS at 
T = 4.2 K) and an important zero-bias background, as in 
Figs.l and 2. These aspects are even more striking if one 
plots the difference between the experimental data and 
the d-wave fit (inset (a) in Fig.l). The shape and magni- 
tude of this difference curve lead to three straightforward 
conclusions. First, the difference curve is local: It is non- 
zero only within the energy range of a few Ao. Beyond 
these energies, there are clearly no important accidents in 
the curves. Second, a simple analysis of this curve shows 
that the number of states found above the cZ-wave fit is 
nearly equal to the number of states found below it, as 
the corresponding areas are compared. Thus, the effect 
is a strong but local modification of the mean-field quasi- 
particle DOS. Third, the general shape of the difference 
curve, in which a local maximum at one energy evolves 
into a local minimum at a higher energy, indicates that 
they both may arise from the derivative of a single peak 
at a characteristic energy situated somewhere between. 

The above arguments, together with the fact that the 
peak-dip-hump appears at the critical temperature, at 
which the pair condensate is formed, indicate that the 
mechanism responsible for these features in the DOS 
should be reflected in the pair function: A (E, <p). An- 
alytically, we choose a Lorentzian centered at ±|-Eo| and 
with the characteristic width 2<5q: 



A(E, tp) = A cos(2<^) 



A S 2 



(E ± E Q ) 2 + S Q 



, (4) 



where A is the maximum value of the non-perturbed 
eZ-wave gap function and A is the normalized amplitude 
(see inset in Fig. 3). The use of (4) for A(E, tp) in (3) gives 
a new quasiparticle DOS having non-trivial characteris- 
tics (Fig. 3). This DOS has the usual , d- wave' singular- 
ity, but at a new energy: Asc, defined by E = A(E, 0). 
In addition, two further extrema arise: a maximum ris- 
ing near Eo — So and one minimum situated near Eo + So, 
both being due to the derivative term. 

At this stage the problem is already reduced to only 
five free parameters: Ao, Eo, So, Ao and the pair- 
breaking r. The fit to the experimental spectra with 
these parameters gives an excellent agreement. Surpris- 
ingly, the best fits do not situate Ao at the quasiparticle 
peak but at some higher energy. Furthermore, the Eo val- 
ues are found to be nearly equal to A for both BSCCO 



and YBCO spectra (solid lines in Fig.l and Fig. 2a respec- 
tively). For simplicity our calculated curves are taken to 
be symmetric with respect to zero energy, whereas the 
experimental spectra are clearly not. However, one can 
improve the quality of the fit considering So and Ao to 
be different for the occupied and empty states in (4) (an 
example is given in Fig. 2b). 

In order to see how far the condition Eo ~ Ao holds 
in general, we a priori set the resonance energy Eo equal 
to Ao in all other fits. In this way, with only four free 
parameters, we succeeded to fit almost all of our data. 
Consulting other reports [fj]-[|Q], we find that our quasi- 
particle DOS should be widely applicable. The fits are 
not as good for the double peaked spectrum in YBCO, 
as in the inset in Fig. 2. The expression (3) allows double 
singularities as shown in Fig. 3d. However, they appear 
slightly smoother than in the TS curves. The possible 
reasons may be that the Lorentzian form (4) is too sim- 
ple to reflect the orthorombicity of this cuprate, and its 
smaller anisotropy. These arguments seem plausible if 
one compares the tunneling data obtained in YBCO us- 
ing STS [|l7| (double peaks) and those from planar junc- 
tions |2(| (simple peaks) . The role of structural or surface 
disorder cannot be excluded either, since in this work we 
observed the spectra of both types in the same thin film 
of YBCO. 

The physical origin of the effect is thus a new cou- 
pling at the energy E ~ Ao perturbing the quasipar- 
ticle spectrum in the superconducting state. In the 
non-superconducting state the interaction is either non- 
existant or inefficient and consequently Ao = 0, and 
possibly So — 0. The gap in the DOS is then 2 A ~ 
2Ao (Fig. 3a or Fig. 3b). One may attribute this non- 
superconducting state to the pseudogap seen at low tem- 
perature J(||l8|]. Indeed, from the latter papers, where 
the pseudogap and superconducting gap were measured 
at the same temperature (4.2 K), one finds the pseudogap 
to be slightly larger than the superconducting gap (inset 
(b) in Fig.l). The values we find for Ao are indeed larger 
than the peak-to-peak gap value A$c extracted from the 
spectra. They are consistent with the low temperature 
pseudogap energy scale. The existence of two 'gap' scales 
(Ao and Asc in our approach) is also suggested in [p7| , 
where Andreev reflection experiments showed the gap of 
28 meV, whereas the cZ-wave fit of the c-axis tunneling 
data gave 19 meV for the gap in the same system. 

In the superconducting state the interaction term be- 
comes non-zero (Ao > 0) and two signatures appear in 
the DOS: An additional peak rises near Ao — <5o and a dip 
near Ao + So- This dynamics comes from the derivative 
term in (3). Qualitatively the phenomenon corresponds 
to shifting states from Ao + So, to Ao — So inside the 
gap Ao (Fig. 3, from (a) to (c) and from (b) to (d)). The 
main quasiparticle peak in the DOS results from the su- 
perposition of two peaks: one at Asc and another at 
Ao — So, and explains why the quasiparticle peaks ap- 
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pear so large when they are not well separated (Fig.l 
and Fig. 2). In this Letter we do not discuss the state 
dynamics in the phase transition to the pseudogap state 
at T = T c . This involves the explicit temperature depen- 
dence of the parameters Aq{T), T(T) and 6o(T), which 
are still controversial. 
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Fig. 3. Illustration of the state dynamics from non super- 
conducting DOS, (a) and (b), to superconducting DOS, (c) 
and (d). In (c) V and Aq are moderate: ZB is small, the peaks 
are single and well pronounced. In (d) V and Ao are large: ZB 
is important, the peaks are broad and doubled, as reported 
in YBCO. For all curves A = 37 meV and So = 17 meV. 
Inset: The shapes of the energy-dependent gap function used. 

As is well known, in strong-coupling approaches to 
conventional superconductivity, the gap is energy depen- 
dent expressing the electron-phonon interaction. This 
perturbs the quasiparticle DOS at characteristic energies 
Tiijph S> Ao and is a second order effect. On the con- 
trary, our gap function is modified significantly near Ao, 
leading to new singularities in the DOS to lowest order. 
In this sense they are an integral part of the supercon- 
ductivity, the relevant parameters being the amplitude 
Aq and the width 28q of the minimum in the gap func- 
tion. The hump has no particular meaning in our model, 
being a consequence of the superposition of the dip with 
the slow decrease of the usual DOS tail. 

In conclusion, we showed that the unusual peak-dip- 
hump features exist in both YBCO and BSCCO, being 
common peculiarities of the quasiparticle DOS in high- 
T c cuprates. Experiments show that these signatures ap- 
pear at T < T c and are intimately related to the super- 
conducting phase. We suggest that the large peaks and 
dips observed in the superconducting state arise together 
from the same interaction, strongly perturbing the usual 
mean-field DOS. This unknown interaction is accounted 



for by introducing the energy-dependent gap function, 
having a minimum at Eq. The numerical analysis of the 
data showed that Eq corresponds to the maximum value 
of the gap Ao for all tunneling spectra in both BSCCO 
and YBCO. The natural consequence of this is the dip- 
hump scaling with the gap peak-to-peak value reported in 
. Finally, the model suggests that the superconduct- 
ing state differs from the 'normal' one by the non vanish- 
ing of dA/dE. These results put strong constraints on a 
possible theory of the superconductivity. 

We are grateful to J. Siejka and Y. Lemaitre for pro- 
viding us the YBCO samples. 
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